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\ Abstract 

, A new nonconforming rectangle element with cubic convergence for the energy norm is intro- 

duced. The degrees of freedom (DOFs) are defined by the twelve values at the three Gauss points 
on each of the four edges. Due to the existence of one linear relation among the above DOFs, it 
turns out the DOFs are eleven. The nonconforming element consists of P2 ffi Spanja; 3 ?/ — xy 3 }. 
We count the corresponding dimension for Dirichlet and Neumann boundary value problems of 
second-order elliptic problems. We also present the optimal error estimates in both broken energy 
and Z/2 (fi) norms. Finally, numerical examples match our theoretical results very well. 

^ . Keywords: nonconforming finite element; optimal error estimates; quadrilateral mesh 

CN 
^> . 

9£ ■ 1 Introduction 

vq 

It has been well known that the standard lowest order conforming elements can produce numerical 
locking and checker-board solutions in the approximation of solid and fluid mechanics problems: see 
for instance [H, 0, 0, 0, EH and the references therein. An efficient approach to deal with this case is 
to employ the nonconformingelement method, which has made a great impact on the development of 
finite element methods [2, 0,1!, 0, 1, EMI MMM,M,M,M,\mMi- 



k> ' To approximate PDEs using a nonconforming element of order A;, one needs to impose the continuity 

of the moments up to order k — 1 of the functions across all the interfaces of neighboring elements. 
This condition is known as the patch test [13] • In two dimensions, the patch test is equivalent to 
the continuity at the k Gauss points located on each interface. This implies that a Pfe-nonconforming 
element, if exists, must be continuous at the k Gauss points on each edge. These points (completed 
with internal points for k > 3) can be used to define local Lagrange degrees of freedom (DOFs) on the 
simplex if k is odd, but this construction is not possible if k is even since there exists a lower-degree 



polynomial vanishing at all the Gauss points |14| . Thus suitable bubble functions are often employed 
to enrich the finite element space. Until now, the triangular nonconforming elements are well studied 
in the literature (see, 13, 14|), but the analysis of their quadrilateral counterparts is less complete. 



Even though the triangular or tctrahcdral meshes arc popular to use, in some cases where the 
geometry of the problem has a quadrilateral nature, one wishes to use quadrilateral or hexahcdral 
meshes with proper elements. For even k, the same trouble exists, that is, there also exists a lower- 
degree polynomial vanishing at all the Gauss points. Again, some bubble functions arc added to the 
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finite element space [2l[. Compared to the triangular case, another trouble for quadrilateral finite 
clement is that the DOFs and corresponding polynomial space do not match. Usually the number of 
DOFs is bigger than the dimension of Pk- For example, for k = 1 and k = 2, the numbers of DOFs are 
4 and 8, respectively, while the corresponding dimensions of Pk are 3 and 6, respectively. Therefore, 
some additional relations must be imposed or some special functions arc added to the finite element 
space such that unisolvency can be satisfied, see 0, El Si]. 



The purpose of this paper is to develop a P3 -nonconforming element on rectangular meshes. We 
define the 12 Gauss points (3 Gauss points on every edge) as the DOFs. To obtain an optimal order 
error estimate, the finite element space must be carefully chosen such that any function in this space 
is a polynomial of degree no greater than 3 on every edge. Meanwhile, we also notice that the values 
on the 12 Gauss points satisfy a linear relation if the degree of a polynomial on every edge is no more 
3, which is a little different from the triangular mesh case. Thus we define our finite clement space as 
P3 © Span{a; 3 y — xy 3 }. Therefore, the number of DOFs is locally 11. We prove unisolvency and define 
three types of local and global bases, one of which is defined associated with vertices and the other two 
of which arc defined associated with edges. Then we derive optimal error estimates for second-order 
elliptic problems in broken energy- and L 2 -norms. Finally, numerical examples are provided, which 
match our theoretical result very well. 

This paper is organized as follows. In Section 2 the P3 nonconforming element is defined on 
rectangular meshes. The dimensions and basis functions for Dirichlet and Neumann problems are 
given. In section 3 interpolation operators are defined and optimal order error estimates are shown. 
In Section 4, numerical results for the elliptic problems are presented. 

2 The P3 nonconforming element on rectangular mesh 

2.1 The ^-Nonconforming Quadrilateral Elements 

Take reference element as square R = [— 1,1] 2 . Denote by Vj,j = 1, ••• ,4, the vertices and by 
53.7-2,33^-1 and g$j the Gauss points on the edges VjVj+i,j = 1,--- ,4, with the identification V\ = V5 
(see Fig. [[]): 

.91 = (-a/375, -l) , .92 = (0, -1), .93 = (^3/5, -l) , 
.94 = (1,-^375) , .95 - (1, 0), g 6 = (l, ^375) , 

97 = (V3/5, l) , 58 = (0,1), 59= (-V37m): 

3io= (-1, a/375) , 311 = (-1,0), g 12 = (-1,V375). 

Denoting by Pk(R) the space of polynomials of degree < k on R, set 

P = P 3 (R) © Span{x 3 y - xy 3 }, with dim(P) = 11. 

The space P will be our nonconforming finite clement space on R with appropriate degrees of freedom 
that will be defined soon. Before proceeding, we notice the following simple result. 

Lemma 1. The following relationship holds: 

3p(-l) + 3p(l) - 5p(v/375) - 5p(— v/375) + M0) =0 e P 3 (M). (1) 



Proof. Let p € Pa(R). Then any fourth order difference quotient of p(x) vanishes. The fourth order 
difference quotient of p(x) at points x% = —1,X2 = — a/3/5, £3 = 0,2:4 = a/3/5, and X5 = 1 can be 
expressed by 

Pjxi) 

~ (Xi - xx) ■ ■ ■ (pa - Xi_i)(xi - x i+l ) ■ ■ ■ (xi - x 5 ) ' 
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Figure 1: The reference rectangle R 



A simple computation derives the desired result. 

As an immediate consequence of Lemma [TJ we have the following proposition. 
Proposition 1. The following relationship holds: 

4(v(32) + <p(gs)) - 5(^(51) + ip(g 3 ) + ip(g 7 ) + tp(go)) 

= 4(^(3.5) + <p(gu)) - 5(^(54) + <p{ge) + <p(gw) + <p(gu)) e 



□ 



(2) 



Proof. Notice that if is a polynomial of degree no greater than 3 on any edge of R. The result follows 
immediately from Lemma [TJ □ 



Lemma 2. Suppose that (p £ P vanishes at the twelve Gauss points g^j = 1, ■ 
R. 



12. Then ip = in 



Proof. Let p £ -P3QR) vanish at vanishes at the twelve Gauss points g 3 r, j = 1, • • ■ , 12. It is easy to 
check that ipi(x,y) — x 2 + y 2 — 8/5 vanishes at the eight Gauss points gsj-2,gsj,j = 1,2,3,4. By a 
simple polynomial division, ip(x,y) can be written as 

ip(x,y) = <pi(x, y){a xy + a\x + a 2 y + a 3 ) + xr\{y) + r 2 (y) 

where r±(y) and r 2 {y) take the following form: 

r i(y) = -2a y 3 + hy 2 + b 2 y + b 3 and r 2 (y) = c Q y 3 + c x y 2 + c 2 y + c 3 . 

All the coefficients aj, bj,Cj are yet to be determined. It follows from <^(±y^3/5, 1) = </Ji(±y3/5, 1) = 
that 

VV^nil) + r 2 (l) = 0, and - + r 2 (l) = 0, 

which implies that ri(l) = r2(l) = 0. Similarly, if follows from (p(±y/3/5, —1) = y>i(±^/3/5, —1) = 
that ri(-l) = r 2 (-l) = 0. Next, it follows from <p(±l, = ipi(±l, V^S) = that ri(±y/3jS) = 

r 2 (iy^3/5) = 0. Since r% and r 2 are polynomials of degree no more 3, which leads to r\{y) = r 2 (y) = 0. 
Hence <p(x, y) must have the following form 

ip(x,y) = (p 1 (x,y)(a 1 x + a 2 y + a 3 ). 

Then the four additional conditions y(53j— 1) = 0,j = 1,2,3,4, lead to (p(x, y) = 0, which completes 
the proof. □ 
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Due to proposition [T] and Lemma [21 we have the unisolvency result. 



Proposition 2. A function tp £ P is uniquely determined by tp(gj),j = 1, 2, • • • , 12, which satisfy the 
relation ([2]) . 

Denote by Mj,j = 1, 2, 3, 4, the four midpoints of four edges. Obviously, Mj is one of three Gauss 
points on jth edge and hence the other two Gauss points on jth edge can also be denoted by and 

MJ . For example, in Fig. [TJ = gi and MJ = g%. For 1 < j < 4, define (pj,(p® + and ip®~ G P 
(see Fig. ^ by 



<Pi (9k) 



(with the identification go = 312), and 



1, = 3.7 - 3,3.7 - 2, 
0, otherwise, 



4, g k = Mf, 

5, gk = Mj, 
0, otherwise, 



and 




<pf (9k 



We then have the following result. 

Lemma 3. Span{^, ipf + , <pf~ ,j = 1,2,3,4} 
span P. 



9k = Mj, 
9k = Mj, 
otherwise. 



Indeed, any eleven of ipj, <pj + , y> ■ ,j = 1, 2, 3, 4, 



The proof of this lemma is similar to that of Lemma 2.3 in [2l|, which will be omitted. 
We are now in a position to state the definition of P3 -nonconforming element on a rectangle as 
follows. 

Definition 1. The P3 -nonconforming element on rectangle is defined by (R, P^, X)_r); where 
• R is a rectangle, 



n 



= P is the finite element space, and 



• ^2 ft = {^(g^jj = 1, 2, • • • , 12, such that Eq. ^ holds for all ip G P^} is the degrees of freedom. 
The following patch test lemma is immediate since any p G P is of degree < 3 on an edge. 

Lemma 4. Let Ej denote the face containing the midpoint Mj and the other two Gauss points Mj 
and MJ. Ifp g P, p{Mj) = p{M+) = p{MJ) = 0, then 



p{x,y)q{x,y) do = 



(3) 



for all q G Qi(Ej 



Remark 1. For actual computation, the local finite element can be alternatively given by {^ ? - =1 (ej >Pj + 



E.L. Ex 



E- E 



Remark 2. In theory, if x 3 y—xy 3 is replaced by x 3 y, xy 3 or x 3 y+xy 3 , the unisolvency like Proposition 
[H holds. But the former two choices are lack in symmetry. As for the third choice x 3 y + xy 3 , it turns 
out to be numerically during the computation of the bases since the corresponding coefficient matrices 
are nearly ill-conditioned (the determinants are near to zero). 
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Figure 2: Three types of basis function, (a) The vertex-based basis function associated with vertex 
V\ . (b) The fist type of the edge- based function associated with edge E\ . (c) The second type of the 
edge-based function associated with edge E\ . 
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Let us proceed to define our P 3 -nonconforming element space. Assume that SI E K 2 is a par- 
allelogram domain with boundary T. Let (^h)h>o be a regular family of triangulation of J! into 
parallelograms Rj,j = 1,2, •• ■ , Nr, where h = maxjjg^ hu with h^, — diam(i?). For each R G 3?h, 
let Fr : R — > R 2 be an invertible affine mapping such that 



For a given triangulation 3?h of fi, let Nv,Ne,Nr, and Nc denote the numbers of vertices, edges, 
rectangles, and Gauss points, respectively. Then set 



In particular, let Ny,N E , and Nq denote the numbers of interior vertices, edges, and Gauss points of 
fieJj, respectively. 

For a function / defined in tt, denote by fj its restriction to Rj, and Ejk the interface between Rj 
and Rk- Similarly, gj^, k = 1, 2, 3, will mean the Gauss points on Tj = dRj n dfl and gjku ^ = 1, 2, 3, 
will be the Gauss points on Ej^ . We are now in a position to define the following nonconforming finite 
element spaces. 

= {ip : — > M.\ip\ji £ P_r, Vi? G 5^, ip is continuous at the Gauss points g G 
^T'g'o = {<? e ./^^(i?) = 0, for all Gauss points gE& h n T}. 

For each vertex Vj £ denote by $h{j) and @h(j) the set of all edges E £ §h with one of the 
endpoints being Vj and the set of Gauss points nearest to Vj among the three Gauss points on E for 
all E £ $h(j)- F° r Mj £ Ej, if g$ and g^ are two other Gauss points and i < k, we also denote these 
two Gauss points by Mj~ and Mj , respectively. We then define the three types of functions in jV c € h ', 
which serve as global bases for the nonconforming finite element spaces. 

Definition 2. The first type of functions are associated with vertices. Define cpj £ jV c S' h ,i = 



R = F R (R), 



and denote 4>r = <fi o F R 1 : R — > K for all £ P, whose collection will be designated by 



P fi = Span-j>,R, G P}. 



%, = {Vi, V 2 , ■ ■ • , V Nv } : the set of all vertices of R £ & h , 
Sh = {Ei,E 2 , • ■ • , E Ne } : the set of all edges of R £ 

= {51:52) ' • ' i9n g } '■ the set of all Gauss points on <§ h £ 3? h 
J( h = {Mi G Ei, i = 1, 2, • • • , iV^} : the set of all Midpoints on S h £ 



1,2,- •• ,N V , by 




1, 9k&&hU) 
0, otherwise 



Next define the second type of functions associated with edges Ej £ $ : define tp^ + £ jY c € h ,j 



1,2,- ■■ ,jV B7 6?/ 




5, = Mj- 
4, ff fc = A/+ 



0, otherwise 



The last type of functions are also associated with edges Ej £ <§ : define ip^ £ jV c £ h ^ j = 1,2, •■ ■ , Ne, 



by 




5, m, Mj 
4, g h = Mj 
0, otherwise 
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Similarly, define the three types of functions which will serve as global basis functions for ^V$ Q 
with those for jV c £ h excluding fj^s which are associated with boundary vertices and tp^ + , tpj ~ 's which 
are associated with boundary edges. 

Now let us present the dimensions for the nonconforming finite element spaces. 

Theorem 1. dim{^^ h ) = N V + 2N E -1- Let(pJ,j = 1,2,--- ,N V and tpf + ,<pf~ ,j = 1,2,--- , N E 
be the functions defined in Definition By omitting any one of these functions, each of the following 
three sets forms global basis functions for JV^ 1 : 

@x = {<pY,<pL--- ,<pYt v -i,<pf + ,<pf-J = l,2,~- ,Ne}, (4) 
^2 ={<pY,j = 1,2,- ■■ ,N v ,cpf\^\... pf-,i= 1,2,- ■■ ,N E }, (5) 



and 



Wj,j = 1,2,--- ,N V , <ff + ,i = 1,2,---N E , <Pi ,<p% , ■ ■ • , <P Nb -i}- ( 6 ) 



Theorem 2. dim(^K<) = +2N%. 28 = {ipY ,i = 1,2,- ■■ ,N^,ipf + ,(pf~,j = 1, 2, • ■ • , N E } 
forms a set of global basis functions for jY^q. 



The proofs of the above theorems are quite similar to those in the literature [21[ , and thus omitted. 
Here we remark that our finite element space is a little different from those in the literature fl8, 21\. 
Those finite element spaces are nothing but conforming element spaces enriched by some suitable 
bubble function spaces. Thus the idea of Fortin and Soulie (3] is not applicable here. 



3 The interpolation operator and convergence analysis 

In this section we will define an interpolation operator and analyze convergence in the case of 
Dirichlet problem. The case of Neumann problem is quite similar and the results will be briefly stated 
with their details being omitted. 

Denote by (•, ■) the L 2 (Q) inner product and (/, v) will be understood as the duality pairing between 
and Hq(£1), which is an extension of the duality paring between L 2 {VL). By || • \\k and | • \k 
we adopt the standard notations for the norm and seminorm for the Sobolcv spaces H k (il). Consider 
then the following Dirichlet problem: 

- V • (aVu) + 0u = /, 0, (7a) 
u = 0, T, (7b) 

with a = (a jk ),atjk,P e L oc (n),j 1 k = 1,2,0 < a*|£| 2 < £*a(a;)£ < a*|£| 2 < oo,£ G M. 2 ,f3(x) > 0,i6 
fi, and / G iJ 1 (Sl). We will assume that the coefficients are sufficiently smooth and that the elliptic 
problem (|7J) has an 7? 4 (f2)-regular solution. The weak problem is then given as usual: find u G Hq{Q) 
such that 

a(u,v) = (/,«), veH^Q), (8) 

where a : Hq(H) x Hq(Q) —> R is the bilinear form defined by a(u,v) = (aVu, Vv) + (f3u,v) for all 
u,v G Hq(Q). Our nonconforming method for Problem states as follows: find Uh G jV c € ] q such 
that 

a h (u h , v h ) = (/, v h ), v h G Jf^Q, (9) 

where 

a h (u,v) = ^ a R( u i v )i 
Resr h 

with an being the restriction of a to R. 
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For a given rectangle R G define the local interpolation operator Hr : W l ' p (R) D Hq(CI) — > 
Pr, p > 1, by 

H-R<l>(9i) = <t>{9%), 

for all Gauss points on the edges of i?. The global interpolation operator : W 1,p (r2)n.ffo(f2) — ► 

is then defined through the local interpolation operator XIr by Hh\n = IIr for all R G 3~h- Since ri/j 

preserves P3 for all R G it follows from the Bramblc-Hilbert Lemma that 

^2 H - ^h(/)\\mn.) +h H-^Wh^r) < Ch k \<t>\ 

R^sr h nesr h (10) 

g w k *{VL) n 2 < k < 4. 

Denote E jk = (dRjDdRk) for all Rj,R k G ST h and by Tj the boundary face (dRjlldD.) of 5^. Then 
define 

A' 1 = {X\X jk = X\ Ejk 6 P 2 (£ ifc ); Aj* + A fcj = 0; A, = A|r, G P 2 (T j )}, 

where P2(E) denotes the set of quadratic polynomials on the face E. Also define the projection 
P h : H 3 / 2 (n) -> A h such that 



v\ aVvj - P h Vj,z) =0 for all z G P 2 (-E) for all E G 

where Uj = and i/J is the transpose of the unit outward normal to Rj. Then we have the following 
standard polynomial approximation result: 

r 1 1/2 

I>J«V, 3 .-P^W ^ C7**" 8/2 IMU. WGW^(Q),A; = 2,3,4,p>3. (11) 

j 

Since — Wfc has zero values at the Gauss points on Ej k for all w G ■yV'&Q and the 3-point Gauss 
quadrature is exact on polynomials of degree no more than 5, the following useful orthogonality holds. 
(See also Lemma [4]) 

Lemma 5. If u G H 3 / 2 (Q), then the following equality holds: 

(P h u. } ,Wj) E]k + (PhUk,Wk) E hj = {PhUj,w. } - w k ) Ejk = for all w G jV^. 
Denote the broken energy norm || • ||^ on 

\\ip\\ h = a h (^,<p) 1/2 for all ip G JT<€ h + H 1 ^). 
We now consider the energy-norm error estimate and first consider the following Strang lemma (26j . 
Lemma 6. Let u G H l (Q) and Uh G ^V'Toq be the solutions of Eq. (jHJ) and Eq. respectively. Then 

( . r \o>h(u, w) — if, if)h , „, 

Assume sufficient regularity such that u G H A . Due to (JTUJ), the first term in the right side of (JT2J) 
is bounded by 

inf ||u-«|| h < \\u-U h u\\ h <ch s \u\ H s+un),l<s<3 (13) 



In order to bound the second term of the right side of (|T2j) which denotes the consistency error, 
integrate by parts elementwise so that 



a h (u,w) - (f,w) = Y "' 1 
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3 

= J2( u J aVu J - PhUj,w - m 3 ) QRAT] 

3 

where rrij G Q2(Rj) is a biquadratic polynomial on Rj. In particular, if rrij is chosen as the Qi 
projection of w in Rj, due to the trace theorem, (fTU)) and (fTTj) . we get 

^(vfaVuj - P h Uj,w- mj) dRj < ^ \Wj a ^ u j - PhUjU^dR^Ww - mj\\ L2{aRj) 

3 3 

^ II vT 3 aVu i ~ p h u i Hl 2 (a^)) ' (Hll w_m J'" 2 



< Ch k -^\\u\\ k ^£ \\w - mjh^Ww - m 3 \\ Hl{Rj} ) 



3\\L 2 (dRj) 
1/2 



(14) 



<Ch k -y 2 \\u\\ k h^\\w\\ h 

= Ch k - 1 \\u\\ k \\w\\ h , k = 2,3,4. 

Thus, we have 

sup { KK^)-(/^)h <^-t h|| fc = 2;3j4 . 

By collecting the above results, we get the following energy-norm error estimate. 

Theorem 3. Let u e H k+1 (Q) n i?o(fi) arl ^ e ^K^o ^ e ^ e solution of (HJ) and ([9]), respectively. 
Then we have 

\\u-u h \\h < ch k \\u\\ H k+i^, k= 1,2,3. 

By the standard Aubin-Nitsche duality argument, the L2(r2)-error estimate can be easily obtained, 
but the corresponding proof is omitted. We state the result in the following theorem. 

Theorem 4. Let u E H k+1 (tt) D Hq(£1) and u h <G be the solution of (J8|) and ^j, respectively. 

Then we have 

\\u -u h \\ < Ch k+1 \\u\\ Hk+ i {n)l k = 1,2, 3. 

Instead of the Dirichlet problem, if the following Neumann problem 

- V • (aVu) + (3u = /, n, (15a) 
v ■ (aVit) + "/u = g, T, (15b) 

is considered, the weak problem ([8]) is then replaced by finding u <E i/ 1 (0) such that 

a n (u,v) = (f,v) + (g,v), v e H^Q), (16) 

where a n is the bilinear form defined by a™ (it, v) = (aVu, Vw) + (j3u, v) + (-fu, v) for all u, v <E i? 1 (fi), 
and (•, •) is the paring between -ff _1 / 2 (r) and iJ 1 ' 2 (r). Thus, the nonconforming method for Problem 
(jTSJ) states as follows: find Uh £ ,jV c i h such that 

al{u h ,v h ) = (f,v h ) + (g,v h ), v h €^ h . (17) 

Then all the arguments given above for Dirichlet case hold analogously Hence one can have the following 
result. 

Theorem 5. Let u <E H k+l (VL) and Uh S ,jV c & h satisfy (fT!)|) and (fTT|) , respectively. Then we have the 
energy-norm and L 2 -norm error estimates: 

\\u - Uh\\h < Ch k \\u\\ H k + i (n) , 

||«-Uh||o < Ch k+1 \\u\\ Hk+Hn) , k = 1,2,3. 
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4 Numerical examples 

In this section we illustrate two numerical examples. First, consider the following Dirichlet problem: 

-Au = /, Q, 

u = o, r, 

where fl = (0, l) 2 . The source term / is calculated from the the exact solution 

u(x, y) = sin(27rx) sin(27ry)(x 3 — y A + x 2 y 3 ). 

Table[T]shows the numerical results, where the error reduction ratios in L2(Q) and broken energy norm 
are optimal. 



h 


DOFs 


||w - Uh\\o 


ratio 


\\u - Uh\\h 


ratio 


1/2 


9 


0.148 




1.759 




1/4 


57 


1.200E-002 


3.62 


0.300 


2.55 


1/8 


273 


4.690E-004 


4.68 


3.051E-002 


3.30 


1/16 


1185 


2.292E-005 


4.35 


3.355E-003 


3.19 


1/32 


4929 


1.279E-006 


4.16 


3.940E-004 


3.09 


1/64 


20097 


7.590E-008 


4.08 


4.78E-005 


3.04 


1/128 


81153 


4.629E-009 


4.04 


5.881E-006 


3.02 



Table 1: The Dirichlet problem: The apparent L 2 and broken energy norm errors and their reduction 
ratios on the quadrilateral meshes. 

Next, turn to the following Neumann problem: 

-Au + u = f, O,, 



with the same domain 17 = (0, l) 2 . The source terms / and g are are generated from the exact solution 

u(x, y) = cos(27rx) cos(27ry)(x 3 — y A + x 2 y 3 ). 

Again, Table [2] shows the numerical results, where the error reduction ratios in L2(f2) and broken 
energy norm are optimal. 



h 


DOFs 


\\u - Uh\\o 


ratio 


\\U - Uh\\h 


ratio 


1/2 


32 


3.850E-002 




0.698 




1/4 


104 


5.217E-003 


2.88 


0.172 


2.02 


1/8 


368 


3.325E-004 


3.97 


2.348E-002 


2.88 


1/16 


1376 


1.917E-005 


4.12 


2.907E-003 


3.01 


1/32 


5312 


1.162E-006 


4.04 


3.616E-004 


3.01 


1/64 


20864 


7.201E-008 


4.01 


4.513E-005 


3.00 


1/128 


82688 


4.491E-009 


4.00 


5.639E-006 


3.00 



Table 2: The Neumann problem: The apparent L 2 - and broken energy norm errors and their reduction 
ratios on the quadrilateral meshes. 
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